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Abstract 

When nonlinearities are "small" there are various ways one can exploit this 
fact — and the fact that the linearized problem can be solved exactly^ — to produce 
useful approximations to the solutions. 

We illustrate two of these techniques here, with examples from phase plane 
analysis: The Poincare-Lindstedt method and the (more flexible) Two Timing 
method. This second method is a particular case of the Multiple Scales approxima- 
tion technique, which is useful whenever the solution of a problem involves effects that 
occur on very different scales. In the particular examples we consider, the different 
scales arise from the basic vibration frequency induced by the linear terms (fast scale) 
and from the (slow) scale over which the small nonlinear effects accumulate. 

The material in these notes is intended to amplify the topics covered in 
section 7.6 and problems 7.6.13-7.6.22 of the book "Nonlinear Dynamics 
and Chaos" by S. Strogatz. 



^Actually, one can also use these ideas when one has a nonlinear problem with known solution, and 
wishes to solve a slightly different one. But we will not talk about this here. 
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1 Poincare-Lindstedt Method (PLM). 

PLM is a technique for calculating periodic solutions. The idea is that, if the linearized 
equations have periodic solutions and < e ^ 1 is a measure of the size of the nonlinear terms 

then: 

I. For any finite time period to < t < to + Tf {Tf > 0), the trajectories for the full 
system will remain pretty close to those of the linearized system (errors no worse than 
0{eTf), typically). 

II. On the other hand, even a small error is enough to destroy periodicity. An orbit that 
"closes on itself" after some time period, will generally fail to do so if slightly perturbed. 
Thus, typically, nonlinearity will destroy most periodic orbits the linearized system might 
have. Some, however, may survive^ — > PLM is designed to pick those up. 

Even if a periodic orbit of the linearized system survives : 

III. The nonlinearity will change (slightly) the shape of the orbit. 

IV. The speed of "travel" along the orbit will be affected by the nonlinearity. In 
particular the period will change (slightly.) 

PLM takes care of these effects as follows; | 

A. The solution is approximated at leading order by the linear solution, but small correc- 
tions at higher orders are introduced to take care of the (small) shape changes. 

B. The linear solution is evaluated at a stretched time, to account for the change in period. 
The two examples that follow illustrate the ideas. 

1.1 Duffing Equation. 

The equation can be written in the form 

X + X + eux^ = , (1.1) 

^That is, if u = u{t) is a periodic solution of the linearized system, then so is au, for any scalar constant 
a. But for only a few values of a will periodicity "survive" the effect of the nonlinearity. 



18.385 MIT (Rosales) Weakly Nonlinear Things: Oscillators. 4 

where < e ^ 1 and u = ±1. This equation is actually a conservative system, with 
(conserved) energy 

E = -x'^ + -x^ + -epx\ (1.2) 
2 2 4 

Thus all orbits for x bounded will be periodic.^ PLM will allow us to calculate corrections 
to the linear period of 27r and sinusoidal orbit shape (for the bounded orbits). 
The PLM expansion is given by: 

x{t) = Xo{T) + exi(T) + e'^X2{T) + • • • , (1.3) 

where Xj = Xj{T) is periodic of period 27r in T and does not depend on e. T = cot is the 
stretched time variable, where 

(J = 1 + ecJi + e^(J2 + • • • , (1-4) 
is a (real, positive) constant to be computed. The nonlinear period is then 27r/ci;. 

Note 1 Xo{T) will be the solution to the linearized problem, so (1-3) will reduce to the 
right answer when e = 0. 

We now proceed as follows: 

• First: Rewrite (1.1) in terms of the new independent variable T, replacing 



/ d 




d . d 


dT 


via 


dt ^ dT 



. Thus: 



cj^x" + x + e = . (1.5) 

• Second: Substitute (1.3) and (1.4) into (1.5) and collect equal powers^ of e. Then require 

that the equation be satisfied at each level in e. Thus we get an equation for each order e^, 

which determine higher and higher orders of approximation in the expansion (1.3), as follows. 

'^Notice that, for v = 1 ALL orbits are periodic. However, for v = —1, orbits where \x\ > are 
not periodic. This follows from looking at the level curves for E in the {x,x) phase plane. Of course, when 
\x\ = 0{e~i), the nonlinear term in equation (1.1) has the same size as the linear terms: the problem is no 
longer "weakly nonlinear". Thus, we should not be surprised if the solution exhibits behavior not close to 
the linearized one. 

^This is the messy part. It means you have to plug (1.3) and (1.4) into (1.5), then do all the products, 
etc. ... so as to end with the equation written as: I {•••} + e {•••} + {•••} + ••• = 0. 
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0(1) equation: 



Xq ~i~ Xq — . 

Clearly then 

Xq = a COS T , 

where a is, at this stage, an arbitrary constant.^ 
0(e) equation: x" + 2u)iXq + Xi + ux^ = 0, 



(1.6) 



(1.7) 



that is: 



x'( + Xi = 2iOiacosT — ua^ cos^ T = 

= l^2u)ia — jua^^ cosT — jua^ cosST . 



(1.8) 



The form of equation (1.8) is typical of all the higher order equations. 

Namely, we get the linear equation for the new term in x at that order — xi 
here — forced by terms involving the lower orders already solved for. 



The solution Xi to (1.8) will be 27r-periodic in T only if the coefficient of the cosT term on 
the right hand side (terms between the brackets) vanishes. This is because this term will 
produce a response in xi proportional to TsinT, which is clearly not periodic. Since we 
are interested in a nontrivial solution (that is a ^ 0) we conclude that: 



Xi = -^ua^ cosST + AcosT + B sinT . 



(1.9) 



where the term marked by the brace in the second equation is the arbitrary homogeneous 
solution, with A and B arbitrary constants. The first equation here determines the first 
frequency correction, in terms of the amplitude^ of the oscillations a, which remains arbitrary 
at this level. ^ We note also that the homogeneous solution in the second equation above 



^In fact, in this case, a will remain arbitrary. There is also a phase shift we could include in (1.7). But 

this is just a matter of where we put the time origin (see appendix A.l). 

^This is typical of nonlinear oscillators: the frequency depends on the amplitude. 

'^That is, no restrictions have been imposed by the expansion on it. In fact, it can be shown that no 

restrictions on a will appear at any level of the expansion. This is because there is in fact a whole one 

parameter set of periodic solutions, which can be parameterized by the amplitude a. 
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amounts to no more than a small change in the amplitude and phase of the leading order 
solution. That is: 

a cos T — > {a + eA) cos T + eB sinT = a cos(r — T) , 
for some a and T. Thus (see appendix A.l) 



Without Loss of Generality: we can set ^4 = = in (1.9). 



(1.10) 



O(e^) equation: 



X2 + 2uiix'( + {2uj2 + Lof) Xq +X2 + "ivx^Xi = 



that is: 



X 



/ \ 9 3 

2 + X2 = (2u)2 + cof) a cosT + —coiua^ cos 3T — —a^ cos^ T cos 3T . 



(1.11) 



16 32 

where cos^ T cos 3T = ^ cos T + ^ cos 3T + ^ cos 5T. Again: X2 will be periodic only if the 
coefficient of the cos T forcing term on the right hand side here vanishes. This yields 

15 



CO2 



-'-2,34 



256 



(1.12) 



and an explicit formula for X2, which we do not display here. Clearly, this process can be 
carried to any desired order (see appendix A. 2). 

In summary, we have found for the solutions^ of the Duffing equation: 



X 


1 Q 

~ acosT + 22^^" cos 


3T+0(e2), 


T 


= cot , 




CO 


~ 1 + ^eua — 25ge a 


+ 0(e3) . 



(1.13) 



1.2 van der Pol equation. 

The equation has the form 

X - ei>{l - x'^)x + X = , (1-14) 

where < e ^ 1 and u = ±1. We use now the same ideas of section 1.1, so that (1.3) and 
(1.4) still apply. Instead of (1.5) we get now 

Lo'^x" + X - euuj{l - x'^)x' = . (1.15) 



^Notice that this is valid only as long as < e ^. When \a\ = 0{e the "corrections" cease to 

be smaller than the leading order and the expansion fails. This agrees with our observations in footnote 3. 
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At 0(1) we get, as before (see appendix A. 3): 

Xq = a cos T. (1-16) 



0(e) equation: x'( + 2coiXq + Xi — z^(l — xI)x'q = 0, that is 



x" + Xi = 2iOiacosT — pa sinT + Ptt^ cos'^ T smT 

(1.17) 

= 2iOia cos T + ua — 1 j sin T + jua^ sin ST . 

To get a periodic solution Xi, both the coefficients of cosT and sinT must vanish on the 
right hand side =^ For a nontrivial solution (a ^ 0) we nnust have^: 

a = 2, (Ji = and Xi = ——ua^ sinST + AcosT + B sinT . (1-18) 

32 ^ V ' 

Note 2 There is an important difference here with the situation in the analog equations 
(1.8) and (1.9). Now both sines and cosines appear on the right hand side of equation (1.17). 
Thus we end up with TWO conditions that must be satisfied if equation (1.17) is to have 
a periodic solution for Xi. These conditions are generally called Solvability Conditions. 
Thus now BOTH a and uji are determined. There is NO FREE PARAMETER left and 
there is just one periodic orbit: the LIMIT CYCLE. 

Since now a is fixed to be a = 2, we can no longer argue that by a slight change in 
the amplitude and phase of Xq, we can set ^4 = = (homogeneous part of the solution, 
marked by the brace above), as we did in (1.10). It is still true, however, that the phase of 
the leading order xq can be changed slightly. We can then use this to conclude (see appendix 
A.3) 



Without Loss of Generality: we can set = in (1.18). (1-19) 



On the other hand, we point out that A remains to be determined. That is, the circular 
part of the limit cycle orbit does not have a radius exactly equal to 2, but rather equal to 
2 + + . . . 



^We could take a = —2 also. This, however, is just a phase change T — > T + tt. Thus, we may as well 
assume a > 0. 
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At the next order (that is, 0(e^)) we will get an equation of the form: 

x'^ + X2 = Forcing. (1.20) 

Again (see note 3) sine and cosine forcing terms on the right will have to be eliminated. 
This will produce two conditions, that will determine both A and i02 uniquely. In X2 an 
homogeneous term of the form acosT will appear,^'' with a and cos determined at O(e^). 
And so on to higher and higher orders. 



A unique solution, up to a phase shift, is produced 
this way to all orders: The LIMIT CYCLE. 



(1.21) 



Note 3 In fact, after some calculation — using (1-16), (1-18) and (1-19) — we can see 
that (1.20) is: 



X2 + X2 = (2u)2 + j^a'^^ acosT + (^ja"^ — 1^ uAsinT 

-^a^ (2 - a^) cos 3T + ^uAa^ sin 3T + cos 5T . 

Thus we conclude 



[1.22) 



1 3 5 

^2 = —7zr;:0'^ , ^ = and X2 = a cos T + -r—zci'^(2 — a^) cos 3T — ^^„^ a^ cos oT , (1.23) 
256 512 3072 

where we recall that a = 2. 



2 Two Timing, Multiple Scales method (TTMS) 
for the van der Pol equation. 

2.1 Calculation of the limit cycle and stability. 

In section 1.1 we basically obtained all the solutions to the Duffing equation (1.1) — 
since we ended up with two free parameters: the amplitude a and an arbitrary phase shift 
r — > r — Tq. On the other hand, in section 1.2 we only obtained the limit cycle solution. 
Now, suppose we want all the solutions to the van der Pol equation (1.14) — this will 

^°With a "f3sinT" homogeneous part of the solution eliminated just as above in (1.19) 
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allow us to determine, in particular, the stability of the limit cycle. The method we introduce 
in this section (TTMS) will allow us to do this. 

The main idea is that, if the solution is not periodic, then we cannot represent it 
with a single solution of the linearized equation (as we did in section 1, with its time 
dependence stretched by co from t to T = cot — to allow for nonlinear corrections to the 
period. For any "short" time period this will be O.K., but over long periods large errors 
may result because they accumulate. To resolve this difficulty we will allow ALL the 
parameters of the linear solution to change SLOWLY in time, so as to track the 
true evolution of the solution. Thus, for equation (1.14), we expand^^: 

X ~ Xo{T,t) + exi{T,t) + e^X2{T,t) H , (2.1) 

where t takes care of the "normal" 27r-periodic dependence induced by the linear solution 
and I r = et I is the slow time (that will allow the linear solution being used to drift (slowly) 
as time evolves, from one linear orbit to the next.^'^) 

Remark 1 Note that now the solution depends explicitly on two times, thus the name 
for the method. In this case the "slow" time is t = et, but in other problems it may be 

T = ('^t — or something else. Figuring out what the exact dependence should be need not be 
trivial and usually requires some thinking: it is related to the rate at which the nonlinearity 
causes drift in the orbits — as opposed to just shape changes. We will talk about this later. 

We now rewrite equation (1.14) in terms of the increased set of "independent" variables 
T and t to obtain (here a dot indicates differentiation with respect to t): 

X + leXj- + c^Xtt + X — evil — x'^)x — e^v{l — x'^)xr = . (2.2) 

Note that the equation is now a P. D. E. ! This method appears to complicate things! How- 
ever, the extra terms are multiplied by e and and so at leading order we only get the linear 
O. D. E. In fact: we will only have to solve linear 0. D. E.'s at each order in the approximation! 

Namely: the orbits in phase space are quite close to the Hnear ones, but the speed at which they are 

tracked is sHghtly different ==> Over long times a big error will accumulate, unless we correct for it. 

^^This is only a first, very simple, implementation. We will introduce a more refined one in section 2.3. 
^^This description, strictly, only applies to xq above. The higher order terms exi . . . are there to account 

for the fact that the nonlinear orbits will have slightly different shapes than the linear ones. 
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As usual, we now substitute the expansion (2.1) into equation (2.2) and collect equal 
powers of e to obtain 



0(1) equation: 



x'^ + Xo = 0. (2.3) 

This is the same as in section 1.2, except that now the arbitrary "constants" in the solution 
of (2.3) will depend on r. We thus have 

xo = Ao{T)e'^ + c.c. , (2.4) 

where c.c. denotes complex conjugate and is complex valued. 

Remark 2 Alternatively, we could write xq = a{T) cost + 6(t) sint, where A = ^{a — ib). 
We cannot now argue, as we did before, that it is O.K. to set b = using the fact that a 
change of time origin t — > t + to is allowed. This is because to has to be constant, while 
setting an arbitrary 6(t) to zero would require to = to{T), at least in principle.^"^ 

Remark 3 The use of complex notation in (2.4) makes life simpler. The kind of expan- 
sions we are doing require at each level of approximation that one expand things like x1 in 
Fourier modes. This is much easier to do with exponentials than with sine and cosines! 



At 0(e) we obtain: 



Xi-\-Xi = —2Xor + T^i^ — Xq)Xo 

= {-2i{^Ao-luAo{l-\Aof))e^'-iiyAy^'}+c.c.. 

This equation is very similar to (1-17), except that now: (i) We are using complex nota- 
tion, (ii) There is no loi term and (iii) A new term in appears because of the allowed 
T dependence. The solution xi will be periodic in t provided the coefficient of the e** forcing 
on the right hand side of (2.5) vanishes. This yields the equation 



dr 



Ao = {l - \Aof) Ao, 



(2.6) 



^^Actually, an argument to set fo = can be made, namely: we expect the solutions of equation (1.14) to 
be basically oscillatory. Thus, they will have maximums and minimums. If we set i = to occur at a local 
maximum, then i; = at i = 0, which yields 6 = 0. But this argument will not work at higher orders. 
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which governs the evolution^^ of the amplitude for the linear (circular) orbits under the 
effect of the weak nonlinearity. 

where a and (p are a real amplitude and phase, respectively, then^^ 



If we let 



-^G2 = and -^a = -i'(4-a'^)a. (2.7) 

These formulas show that the orbits in the phase plane are nearly circular, with a slowly changing 
radius a that evolves following the second equation in (2.7) and a limit cycle for a = 2. In 

particular: 



[For 1^ = 1 the limit cycle is stable and it is unstable for u = — l" (2.8) 



If we let = eu in (1.14) and write the equation as 

X- ii{l-x'^)x + x = 0, (2.9) 

then we see that our calculations here show that at = we have a bifurcation, with an 
exchange of stability between the limit cycle and the critical point at the origin. 



// < 0. Unstable limit cycle and stable spiral point. 

> 0. Stable limit cycle and unstable spiral point. 
[1 = 0. Center with continuoum of periodic orbits. (There is no limit cycle.) 



(2.10) 



2.2 Higher orders and limitations of TTMS. 

We us now finish the 0(e) calculation and solve equation (2.5) using (2.6). We have 

X, = i^^ipAy-'' + Ai(T)e^*| + c.c. , (2.11) 

where Ai is complex valued. 

Let us now continue the expansion to one more order, as there is an important detail to 
be learned from doing this. 



Drift in phase space 

Since this shows that (p is a constant, we could have taken fo = in remark 2 ! 
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The O(e^) equation is: 



X2 + X2 = —2Xit- — Xorr + + UXor — UX^Xi — 2vXqXiXq — VX^Xor 

= {-2i - ^i^A, + i^\Al\A^ + ^uAlAl + I'luA', - \iA!^ (2.12) 
-\iv {\Al\ + e** + (. . .)e-^ + (. . ^e'^] + c.c. , 



where 



and ^4* denotes the complex conjugate of A^. Thus, to avoid secular terms 
in X2 (namely: terms proportional to te**, that destroy the periodicity in t) the coefficient 
of e'* on the right hand side of this last equation must vanish. Thus 

A[ - Ka, + V \Al\ A, + \vAlA\ = -UvA, + Uj^', + Uv A,) - A, . (2.13) 

This is a rather messy equation. We do not aim to solve it here; but only to analyze its behavior 
for T large. 

I Assume = 1;^ In this case the limit cycle is stable and, for r large — see equation 
(2.7) — ^0 ^ e*'^, for some constant Then equation (2.13) reduces to 



This is much simpler and can be solved explicitly 



17 



A^ = Cie"^ + iQ 



tip 



(2.14) 



(2.15) 



where Ci and C2 are real constants. This means that the solution of equation (2.13) 
will behave, for large r, like 



16 



-ITS 



lip 



(2.16) 



This is "bad". Notice that the expansion (2.1) for the solution of (1.14) — use equations 
(2.4) and (2.11) — is 

X - 2 Re (Aa{T)e'^) - Im [Al{T)e^'*) + 2eRe (Ai(T)e^*) + • • • . 

But, when er = 0(1) the second term in the expansion will not be small at all (as 

eAi - -^iere'^) ! Thus 



The two timing expansion (2.1) is only valid as long as |t| <C e ^ (2-17) 



I'^Write = ze''^. Then z' + Re{z) 



18.385 MIT (Rosales) Weakly Nonlinear Things: Oscillators. 13 

This is pretty typical for TTMS expansions: Usually they are valid for a time range 
where the "slow" time can be taken large — but not arbitrary large. Beyond some e~^, for 
some p, they fail. 

In the current situation (2.17) is not terribly upsetting. It still allows us to take r fairly 
large. Once r is large and the limit cycle is reached =^ can switch to the expansion in section 
1.2 !! 

However: suppose that (2.17) makes us terribly unhappy, for whatever reasons. Then 



Can we fix the problem posed by (2.17)? (218) 



The answer to this question is YES, but first we must understand why (2.17) 
occurs! This is clarified next; for simplicity we CONSIDER ONLY the STABLE 
LIMIT CYCLE case, when u = 1. 

Note 4 Equations (2.1)-(2.7) lead to an approximation of the limit cycle (for large t, so 
that Aq ~ e*'^ ) given by 

x - 2 i?e(e'(*+'^)) = 2 cos(t + cp) . (2.19) 
On the other hand, the PLM calculation of section 1.2 tells us that we should use 

X - 2 cos(wt + ip) = 2 i?e(e'('"*+'^)) , 

where io = 1 — Y^e^ + • • •. Now, since (expand in Taylor series) 

gi{u:t+v) ^ gi(t+v)g-i^e^t+... ^ ^iit+v) _ Li^^^iit+V) + . . . ^ ^2.20) 

16 

we see that the error in (2.19) is — ^ze^te*^*"'"'^^ + • • • , which is precisely the "bad" behavior 
arising in Ai earlier in equation (2.16). Thus 



The TTMS expansion goes bad because it does not properly take into account 
the fact that the nonlinearity affects the phase — i.e. the position along the 
linear orbit of the solution. 



(2.21) 



• It follows that, to achieve (2.18) we must fix the problem pointed out by (2.21). 
THIS WE DO NEXT. 
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2.3 Generalization of TTMS to extend the range of validity. 



Let (/) be the phase of the solution — namely: its position along the orbit — and 
its angular velocity. The phase increases with time and, for the linearized equation, we 
have 

^(f) = uj = l. (2.22) 

However, once nonlinear effects kick in, there is no reason for co to remain equal to 1, 
or in fact even constant! 

Now, when considering a periodic orbit, as long as u) is approximated by its correct 
average value things will be O.K. (as then errors will not accumulate over time). This is 

J with I (J = 1 + eijQi + • 



what PLM does, by taking 



iOt\ 



. We cannot use this 



idea of PLM in TTMS, because now the orbit (thus the average value of w) varies slowly 
as time changes. We must then allow cj to be a function of r. Thus 
To fix the type of problem discussed in the previous section 2.2 we must replace 
the expansion (2.1) by a subtler type, where the phase (fast time) itself is to be determined. 
Generally we must deal then with expansions of the form 

x-Xo(r,0) + fXi(r,0) + e2X2(r,0) + -- - , (2.23) 

where 27r-periodic dependence on the phase (j) is required, t = et and 

^(j) = uj = 1 + ea)i(T) + 0)2(1") + • • • . 



This amounts to writing: = ^ (r + e (pi(r) + 02 (t) + ■■■), where 

When no r dependence is allowed, this reduces to PLM. We will not carry out the details of 
this expansion here — they are quite messy and some technicalities are involved in selecting 
the (Jj's so that the Xj's behave "properly" as functions of r (that is, no secw/ar growth in r 
occurs). On the other hand, in the particular case of the van der Pol equation (1.14), when 
the limit cycle is stable^^: all solutions eventually approach the limit cycle, and they do so 
on time scales where r <^ (as follows from our results in section 1.2). Thus, as long as 
no cumulative errors occur in tracking the limit cycle, there should be no problems. We can 
conclude thus, without doing any calculations, that: 



i^That is, v=l. 
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For equation (1-14), in the case u = 1: 

• The Lo/s in equation (2.23) are constant and equal to the values calculated for the 
expansion in section 1.2. 

• The functional form of Xo{t, (/)) in equation (2.23) is the same as that we obtained 

for Xq{tJ) in equation (2.1), with t replaced by 0. That is; A'o(r, 0) = .i"o(r, 0). 

In particular, note that from this we learn that the TTMS approximation for the behavior 
of the van der Pol equation is quite good. The secular growth displayed by Ai in equation 
(2.16) for very long times is nothing to worry about. It is simply a manifestation of the 
fact that we have some small (very small, 0{e^ j) errors on the velocity at which the solution 
moves along the limit cycle, but of nothing else. No important qualitative or quantitative 
effect is missing. 

Note 5 Other ways to fix the problem in (2.17) can be devised. For example, some people 
advocate introducing ever slower time scales, such as eH, e^t and so on — in addition to 
the et of equation (2.1). This is not a good idea, unless the problem truly depends 
on that many scales! For example: if the difficulty arises because the true slow time 
dependence^^ is on something like (say) j:^t and not et, then this "lots of scales" approach 
will just complicate things for no real gain at all. For an expansion to be useful, it has to 
zero into the real behavior of the solution. The aim of doing an asymptotic expansion should 
be to learn something useful about the solution, not to produce a massive amount of algebra 
(even if this is, sometimes, an unfortunate byproduct, it is not the aim). In particular, 
producing an "approximation" that fools us into believing that the solution depends on very 
many different time scales (when in fact it does not), is exactly opposite to this objective. 



^^Notice that the van der Pol equation is exactly an example of this type. 
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A Appendix. 
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A.l Some details regarding section 1.1. 

Generally, asymptotic expansions — like the ones in these notes — require at each level 
the solution of a linear equation with some forcing made up from the prior terms. The 
solution of this linear equation is then required to satisfy some condition (periodicity in the 
examples here) and this imposes restrictions on the forcing terms. These restrictions are 
then used to determine free parameters, slow time evolutions, etc. 

When solving the linear equations in the expansion, it is very important to include in the 
solution ALL the free parameters consistent with the conditions imposed on the solution. 
This is because parameters that are "arbitrary" at some level, may later be needed to satisfy 
the restrictions at a higher order. ^° Failure to include a particular parameter — which boils 
down to setting it to some arbitrary fixed value — will typically cause trouble at higher 
order, when a restriction on a forcing term will be found impossible to satisfy. 

On the other hand, practical considerations dictate that we carry as few free parameters 
in a calculation as feasible. Thus, one must always look at the equations involved and ask 
if there is some argument that would allow for the elimination of a parameter — but never 
must one eliminate a parameter without a good reason. 

Consider now equation (1.1) — or (1.5). This equation is invariant under time trans- 
lation: if a: = X{t) is a solution, then so it is a: = X(t — to). Thus, we can always pick the 
origin of the time coordinate to simplify the solution and eliminate parameters. 
For example: The general solution of (1.6) is: acos(r — Tq), where a and Tq are constants. 
But the invariance under time translation shows that we can set Tq = 0. 
Furthermore: At the level of (1.9) we know that in fact a is arbitrary. Then, since A and 
B in (1.9) amount to making small 0(e) changes to a and Tq at the 0(1) level — thus they 
are not true "new" free parameters — we can again set ^4 = = 0, as in (1.10), without 
any fear. 

^°For example, in section 1.2, the amplitude a in (1.16) is eventually set to a = 2 in (1.18). 
Conversely: if an expansion fails at some level, one should always check to see if somehow an important 

degree of freedom (some parameter) was ignored! 
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At any level 0(e") in the expansion, for n > 1, we can 
take Xn in (1.3) with NO cosT or sinT components. 



(A.l) 



A. 2 More details regarding section 1.1. 

It is clear that, in the expansion of section 1.1, the 0(e") equations — for n > 1 — have the 
form 

n 

-\- Xn = Pn{Xo, ■ ■ ■ , Xn-l) — ^ a£Xj^_^ , (^-2) 

t=l 

where P„ is a cubic polynomial and the a^s are constants defined by = X]^o Thus 
Qfo = 1, Qfi = 2iOi, 0L2 = 2(^2 + Wi, aa = ScJa + 2ci;iCi;2, ■ ■ ■ ■ In general we can see that 
a„ = 2^',i + /fi(-^'i- ■ ■ ■ -■^Vi-i)- where /„ is a quadratic polynomial. 



Because Xq is even, the forcing on the right hand side of (1.8) is also even. Then (1.10) 
gives Xi even. The same type of argument shows then that X2 is also even. More generally, 
one can show using (A.l) that all the a;„'s are even. 

Now, the condition on (A. 2) to get Xn periodic in T is that the right hand side should 
not have any forcing proportional to either sinT or cosT. But the right hand side is even, 
thus there is NO sinT forcing ever. On the other hand, the coefficient of the cosT forcing 
has the form: 2au)n + Gn{a, uJi, . . . , a;„_i), where G„ is some polynomial function. Thus, one 
can always choose a;„ so as to make the coefficient of cosT vanish. We have thus 
shown that 



The expaiision in equation (1.3) works up to any order. (A. 3) 



A. 3 Some details regarding section 1.2. 

Equation (1.14) is invariant under time translation. Thus, just as we did in appendix A.l, 
we have a phase to play with and can use to eliminate parameters. 

We used this fact in (1.16) to eliminate the sine component in Xo{T). But now a is no 
longer a free parameter in the solution, as equation (1.18) shows that a = 2. Thus, in order 
to eliminate spurious parameters in Xi{T) (from the two - A and B - that appear in (1.18)), 
we only have a phase to play with. 
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Since 2cos(r— ^eB) = 2 cosT + eB sinT + . . ., it follows that a small phase change can be 
used to eliminate B in Xi{T) as given in (1.18). But A cannot and should not be eliminated 
from the formula. In fact, at O(e^) the solvability requirement on the equations (periodicity 
of a:2(r)) will determine A in the same fashion that a = 2 followed from the 0(e) equation. 
At this level it will be possible to argue that no term in sinT is needed in X2{T), but a 
term acosT must be kept (with a determined at 0{e^)). Clearly the same pattern will be 
repeated over and over. In this fashion the expansion can be continued to any desired 
order. 



